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ABSTRACT 

A fully-asynchronous network with one target sensor and a 
few anchors (nodes with known locations) is considered. Lo- 
calization and synchronization are traditionally treated as two 
separate problems. In this paper, localization and synchro- 
nization is studied under a unified framework. We present 
a new model in which time-stamps obtained either via two- 
way communication between the nodes or with a broadcast 
based protocol can be used in a simple estimator based on 
least-squares (LS) to jointly estimate the position of the tar- 
get node as well as all the unknown clock-skews and clock- 
offsets. The Cramer-Rao lower bound (CRLB) is derived for 
the considered problem and is used as a benchmark to analyze 
the performance of the proposed estimator. 

Index Terms — Clock synchronization, clock-skew, clock- 
offset, localization, wireless sensor networks. 

1. INTRODUCTION 

Localization and clock synchronization are two key compo- 
nents of any self-organizing location-aware wireless sensor 
network (WSN). A WSN enables distributed information pro- 
cessing tasks like data sampling, information fusion, and other 
time-based tasks Q). Every node in the network has an au- 
tonomous clock. These individual clocks in a WSN drift from 
each other due to imperfections in the oscillator, aging and 
other environmental variations. It is essential to calibrate these 
imperfections from time to time to achieve a network-wide 
time coherence. A plethora of algorithms for clock synchro- 
nization can be found in |2j. For the data to be meaning- 
ful, the location where the data is acquired is often required. 
Computing the location of the nodes is commonly called lo- 
calization, and is a well-studied topic 0. 

Even though localization and clock synchronization are 
tightly coupled, traditionally they are treated as two sepa- 
rate problems. Recently, for an anchorless and a fully asyn- 
chronous network, a global least-squares (GLS) estimator based 
on a two-way time-stamp exchange protocol for joint clock 
synchronization and ranging has been proposed in 0. Ex- 
ploiting the broadcast nature of the wireless medium, an asym- 



metrical time-stamping and passive listening (ATPL) protocol 
was proposed in [5 1 for joint clock synchronization and rang- 
ing. Subsequently, the estimated pairwise distances can be 
used as an input to the least-squares (LS) based range-squared 
method for localization. Joint estimation of the position and 
the clock parameters of a sensor based on the two-way time- 
stamp exchange protocol has been considered in [6|, where 
an synchronous network is considered where only the sensor 
node suffers from clock-skews and clock-offsets and the an- 
chors are assumed to be synchronized. In Q, localization 
of the sensor node in a fully-asynchronous network has been 
studied, where the main focus is on localization but also cer- 
tain approximations of the clock parameters are required to 
solve the problem. 

In this paper, we again consider fully-asynchronous net- 
work consisting of one sensor and a few anchors, and inves- 
tigate localization and clock synchronization under a unified 
framework. We propose a linear data model for joint localiza- 
tion and clock synchronization. The data model is generic in 
the sense that it can be used with either the two-way ranging 
protocol or the ATPL protocol. This is used in an estima- 
tor based on LS to jointly estimate the position of the sensor 
node and all the unknown clock-skews and clock-offsets. The 
Cramer-Rao Lower Bound (CRLB) is derived for the consid- 
ered problem and is used as a benchmark to analyze the per- 
formance of the proposed estimators. 
Notation: 

Upper (lower) bold face letters are used for matrices (column 
vectors); (-) T denotes transposition; (0) refers to element- 
wise matrix or vector product (division); (.) 02 denotes the 
element-wise matrix or vector squaring; bdiag(.) a block di- 
agonal matrix with the matrices in its argument on the main 
diagonal; ljy (Ojv) denotes the iV x 1 vector of ones (zeros); 
Ijv is an identity matrix of size N; E(.) denotes the expecta- 
tion operation; Cg> is the Kronecker product. 



2. NETWORK MODEL 

We consider a fully-asynchronous network with M anchors 
and one sensor (node 0) as shown in Fig.[T] We assume one 
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Fig. 1: An illustration of the network model, together with 
the know and unknown parameters. Light shaded lines refer 
to the passive listening links. 

of the nodes has a relatively stable clock oscillator and is used 
as a clock reference. All the other nodes suffer from clock- 
skews and clock-offsets. The network model considered here 
is the same as the model considered in 0. 

All the nodes are distributed over an Z-dimensional space, 
with I — 2 (plane) or I — 3 (3-D space). Let the vector 
Xj G M. lxl denote the coordinates of the ith node. The co- 
ordinates of all the anchors are collected in a matrix X = 
[xi, X2, . . . , Xm] € M ixM . The unknown coordinates of the 
sensor are collected in Xq. The distance between the ith node 
and the jth node is denoted by 



Xj -Xj 2 = 



I 2 - 2xF> 



(1) 

Let U be the local time at the ith node and t be the refer- 
ence time. We then assume that the relation between the local 
time and the reference time can be given by a first-order affine 
clock model [4], 



t = OL l t l + ^ 



(2) 



where cj s ; 6 R + is the clock-skew, fa £ K is the clock-offset, 
cti = uj^ 1 and Pi = —uj~ 1 (f)i are the synchronization parame- 
ters of the ith node. Without loss of generality, we use anchor 
M as absolute time reference, i.e., [ljm,<I>m] — [LO]. The 
unknown synchronization parameters are collected in a = 
[a ,ai, . . . ,a M -i] T and/3 = [p , Pi, ■ ■ ■ , Pm-i] T ■ The un- 
known clock-skews and clock-offsets are respectively given 
by 

ui = 1m oc and <fi — —f3 a. (3) 

Nodes in the network can communicate with each other 
either via a two-way communication protocol [4 1 or an ATPL 
protocol ||5] as illustrated in Fig.Q] 



3. PROBLEM FORMULATION 

In this paper, we focus on the two-way communication proto- 
col for deriving the data model. The transmission and recep- 
tion time-stamps are recorded independently at local time co- 
ordinates both during the forward and the reverse links. The 
time-stamp recorded at the ith node when the kth iteration 
message departs to the jth node is denoted by and upon 
arrival of the corresponding message, the jth node records the 

(k) 

time-stamp ' . For the sake of generality, we do not put any 
constraints on the sequence of forward links and reverse links 
or the number of time-stamps recorded |@]|5]. 

The time-of-flight for a line-of-sight (LOS) transmission 
between the ith node and the jth node can be defined as = 
v^dij, where v denotes the speed of a wave in a medium. 
Using (0, Tij can be written in terms of the local clock coor- 
dinates as 



(a J T^+P 3 )-(a t T^+p i ) + 



(fc) 



(4) 



(k) 

where ' denotes the aggregate measurement error on the 
time-stamps. 

In all there are K time-stamps recorded at each node and 
the time-stamps recorded at the ith node are collected in t ^ = 

[T-P,T?P, T, ( f } 1 T e M A ' xl . The direction (forward or 

(k) (k) 

reverse) of the kth link is denoted by e\j , where = 1 for 

(k) 

transmission from node i to node j and e\-' = — 1 for trans- 
mission from node j to node i. The direction information is 

J K )]T c ]u>Kxl 
ij > "13 ) ■ ■ • > e ij \ £ K 

(1) (2) (K) ]T 



collected in a vector ey = [ey ,e. 
The error vector is denoted by ny = [n 

TaKxl 



For the sake of exposition, we consider the example of a 
network with M = 3 anchors and one sensor (node 0) with 
a of two-way communication protocol between each of the 
sensor-anchor pairs. Let the clock parameters corresponding 
to the ith node be collected in a vector Cj = [a,-,/3j] and 
To = [to,i,t 0) 2,...,to,m] T e K Mxl . The pairwise dis- 
tances of the sensor to each anchor will be do = vtq. We 
can now write (0]l for all the K time-stamps collected in a 
matrix-vector form shown in (|5) on top of this page. Mov- 
ing the known columns corresponding to C3 = [1, 0] T (clock 
reference) to one side, we can re-write (|5} shown in (]6} on 
top of this page. The generalization of (O for any M > 2 is 



straightforward. 

In case we adopt the broadcast based ATPL protocol, the 
matrix A will have additional rows corresponding to the pas- 
sive listening links 0. The detailed derivation of the linear 
data model for the ATPL protocol can be found in [5 1. 

The generalized linear model for either the two-way com- 
munication or ATPL protocol can be written as 

A0 = t + n (7) 
where A e rkmxsm^ q g R 3A/xi^ t £ rkmxi^ and n £ 
R™ xl with the structures detailed in © for the two-way 
communication protocol or in [5 1 for the ATPL protocol. 

The aim of this work is to estimate the position Xo of 
the target node along with all the unknown clock parameters. 
The position of the target node xo can be computed using the 
range estimates obtained by solving (|7). This is presented as 
a two-step approach in the next section, where in the first step 
we estimate the unknown clock parameters and the pairwise 
distances of the sensor to each anchor, and use this estimated 
pairwise distances to compute the position of the target node 
in the second step. Alternatively, we can formulate a single 
estimation problem to jointly compute the position of the tar- 
get node as well as all the unknown clock parameters and this 
is the main contribution of this paper. 

4. TWO-STEP APPROACH 

In the two-step approach, we first solve for all the unknown 
clock parameters and the pairwise distances of the sensor to 
each anchor and then use the range estimates in a LS estimator 
to compute the location. 

4.1. Joint synchronization and ranging: step I 

For K > 3, matrix A is tall and left-invertible. Hence, the 
unknown parameters in can be estimated using LS, i.e., 

= (A T A)" 1 A T t. (8) 
Subsequently, the clock-skews u>, clock-offsets <j> can be ob- 
tained using the relation in ||3), and the pairwise distances of 
the sensor to each anchor using the relation do = vtq- 

4.2. Localization from estimated pairwise distances: step II 

Pairwise distances form a major input to any localization scheme. 
Using the pairwise distance estimates obtained in ((8), the co- 
ordinates of the sensor node can be estimated using range- 
squared localization algorithms. 



Let us define a vector q 



|xi|| 2 ,||x 2 || 2 , 



|xm|| 2 F G 



l Mxl . Using ([TJ, we can write the pairwise distance of the 
sensor to each anchor in a vector form as 



d 0d o 



q- 

Xp 



2X T x 



|xo|| 2 l 



M 



(9) 



wherep = [xg\ ||x |l 2 ] T andX = [-2X T , 1 M ] G R Mx( - l+1 l 
Subsequently, the coordinates of the sensor can be estimated 



using LS as follows 

P = (X T X)- 1 X T (d o 0d o -q) 



(10) 



provided M > I + 1 such that the matrix X is tall. The 
anchor positions can be designed such that the matrix X is 
left-invertible. 

5. THE JOINT ESTIMATOR 

Although clock synchronization and localization problems are 
tightly coupled, they have a non-linear relation as can be seen 
in (|9]l. In case we want to localize the sensor and also estimate 
all the unknown clock parameters in one linear step, we have 
to linearize the relation between the clock parameters and the 
position. 

In order to do such a joint localization and synchroniza- 
tion, we exploit the linear relation between the range-squared 
measurements and the coordinates of the sensor. Instead of 
squaring the estimated pairwise distances in the second step 
as in Section 14.21 we linearize the problem by squaring the 
linear model in ©, and this is the main contribution of this 
paper. A unified framework for localization and synchroniza- 
tion is essential for applications in WSNs such as joint track- 
ing of the position and the clock-parameters, for e.g., using a 
standard Kalman filter. 

Hadamard squaring the data model in (|7) would result in 
a linear model and is given by 

(A0) (A0) = (t + n) (t + n) (11) 
and can be further simplified to 

(A T o A T ) T (0 0) = (t + n)0(t + n) (12) 

The linear model in (fT2l does not have a unique solution as 
the system matrix (A T o A T ) T is a KM x 9M 2 matrix which 
is generally fat and thus not left-invertible. However, this can 
be solved with certain approximations of the clock parameters 
as in J7J- 

An alternative approach to linearizing the problem is by 
taking the Kroneckor product of the measurements, i.e., 



(A0) (A0) = (t + n) (t + n). 



(13) 



Using the matrix property PC QE = (P Q)(C E), 
we can further simplify ([T3~t to the following linear model 



a e 

(A A) (0 0) = (t + n) (t + n) 



(14) 



1 1 + n ( 



n+t0n+n0t 

where A e rK 2 m 2 x9M 2 ^ and w g r k 2 m 2 xi is the new er _ 
ror vector. If the matrix A is full column-rank, then it follows 
that the matrix A is also full column-rank. 



We now introduce two new variables to resolve the clock 
parameters without ambiguity after the squaring operation. 
For the ith node, we define the variables 



7* 



and Si = ouPi 



(15) 



and collect parameters corresponding to the ith node in the 
vector c, = [7i,<5i] T . The clock-skew uj; L 6 M. + is always 
positive and the clock-offset G M can be either positive or 
negative. As a result, recovering clock-offsets from f3f with- 
out ambiguities would be difficult. Hence, we make use of the 
cross-term Si = ctifii to recover the clock-offset. For all the 
nodes in the network we have c = [cq , cf , c^, . . . , cJ f _ 1 ] T . 
Let us define a permutation matrix II G R 9m2 x9A/2 that sorts 
the entries of 6, such that 110 = [c T , Tq 2T , z T ] T . Here, the 
entries of the vector z G R LzXl with L z = 9M 2 - 3M, con- 
sist of the nuisance parameters excluding c and rf 2 from 9 
and is of less interest. 

We can now re-write (fl4t as follows 



An T (s 5 c 



- 2 s d d® 2 + s z z) = t + w 



(16) 



where S e , S d , and S z are the selection matrices to select 
columns of AII T corresponding to c, dj^ 2 , and z, respec- 
tively. Substituting (O in ([ToD . we get 



An T (Scc + ^ 2 s d x aP + s zZ ) 



t - An J z^ 2 S d q + w. 

(17) 

We next collect the unknowns in the vector ijj = [c T , p T , z T ] T 
G R Lxl where L = 2AI + I + 1 + L z and the columns cor- 
responding to the unknowns in the matrix 

A = [An r Se,^ 2 An T s d x, An T s z ] e r k2m2xL , 

- - —O X T ^ 18) 

and the measurements in the vector t=t — v All S d q G 

^K 2 M 2 xl 

The generalized linear model for joint localization and 
synchronization is then given by 

A-0 = t + w (19) 
The unknown parameters in tp can be estimated using LS, i.e., 

# = (A T A)- 1 A T t. (20) 

Hence, the unknown position xo is obtained by solving < f20b 
and the unknown clock-skews and clock-offsets can be ob- 
tained using (Tl3T > and (O without any ambiguities. 

Alternatively, a weighted least-squares (WLS) estimator 
instead of ( TTOb taking the estimation error in (|8) or a WLS es- 
timator instead of (l20b pre-whitening the noise w is possible. 
However, this is not further detailed in this paper. 



unbiased estimator tp it follows from the CRLB theorem that 
K(xpxp T ) > F _1 where F is the Fisher information matrix. 
If the error vector n is Gaussian distributed with a variance 
a 2 , then F can be computed as F = o~ 2 J T J, where J is the 
Jacobian matrix given by 



9(A0-t) 



[Ju 



10 



*KMx(2M+l) 



dip 

with sub-blocks 

= -A(S Q -S/3 1km4> T ) 

= -AS/3 1kmv t , 
J xo - i/^TS^D, 



(1kmu> ) 



T\Q2 



(21) 



(22) 



where S a , Sp, and St- are selection matrices to select the 
columns of A corresponding to a, (3, and To, respectively. 
The M x I derivative matrix D is defined as 



[D" 



• ■j 



ado 

<9x 



l|xo-Xi|| 2 



(23) 



7. NUMERICAL EXAMPLE 

A network with one target sensor and 5 anchors is considered. 
Both the target node and anchor nodes are deployed randomly 
within a range of 100m. The clock-skews u> and clock-offsets 
<f) are uniformly distributed in the range [1 — lOOppm, 1 + 
lOOppm] and [—Is, Is], respectively. We use an observation 
interval of 100 s during which the clock parameters are as- 
sumed to be fixed. We use v — 300 m/s and record K = 10 
time-stamps. The error vector is assumed to be Gaussian dis- 
tributed with a variance a 2 . The simulations are averaged 
over 1000 independent Monte Carlo experiments. 
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Fig. 2: RMSE of the estimated sensor coordinates. 



6. CRAMER-RAO LOWER BOUND 

We now derive the CRLB for jointly estimating the clock- 
skews ui, the clock-offsets <j>, and the coordinates of the sen- 
sor node xo, i.e., xp — [o; T , cp T , x.q] t based on 0. For an 



In this paper, we analyze the performance of the proposed 
estimators in terms of the root mean square error (RMSE) of 
the estimated sensor position, clock-skews, and clock-offsets 
for different values of a 2 . We provide the results for a) two- 
way communication protocol [4 | b) ATPL protocol 0. 




Fig. 3: RMSE of the estimated clock-skews. 
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Fig. 4: RMSE of the estimated clock-offsets. 

Fig. |2] shows the RMSE of the estimated sensor position 
computed using a) the two-step approach, i.e., LS estimator 
( fTol ) with the range estimates of (HJ as the input, and b) the 
proposed joint estimator. The root CRLB (RCRLB) is also 
provided for both the considered protocols. The ATPL pro- 
tocol performs better than the two-way communication pro- 
tocol due to the additional passive listening links |5|. How- 
ever, both the estimators for localization are inaccurate as the 
dependencies with the nuisance parameters are not consid- 
ered, and this can be resolved using a constrained WLS solu- 
tions Q. 

Fig. [3]and Fig. |4]show the RMSE of the estimated clock- 
skews and clock-offsets. The ATPL protocol again performs 
better than the two-way communication protocol. In addi- 
tion, both the estimators for clock-skews and clock-offsets are 
asymptotically efficient, and meet the CRLB. 

The location of the target node can be obtained with a two- 
step approach using the range estimates obtained from joint 
synchronization and ranging. Alternatively, we can formulate 
localization and synchronization under a unified framework 
as a single linear problem. However, this results in a larger 
system to solve and is computationally less attractive than the 
two-step approach. The linear model of the joint synchroniza- 
tion and localization can be used for joint tracking of the clock 



parameters and the position which is an important application 
in a WSN. 

8. CONCLUSIONS 

We have considered a fully-asynchronous network with one 
sensor and a few anchors. In this paper, we have addressed 
a problem in which we estimate all the unknown clock pa- 
rameters as well as the position of the target node. Loca- 
tion of the node can be estimated with a two-step approach 
using the range estimates. To avoid this two-step approach, 
we have proposed a generic linear data model for joint lo- 
calization and clock synchronization. An estimator based on 
LS to jointly estimate the position of the target sensor node, 
along with all the unknown clock-skews and clock-offsets has 
been presented. The proposed estimator for clock-skews and 
clock-offsets is asymptotically efficient and meets the CRLB, 
however, the position estimates do not asymptotically achieve 
the CRLB. 

9. REFERENCES 

[1] N.M. Freris, H. Kowshik, and P. R. Kumar, "Funda- 
mentals of large sensor networks: Connectivity, capacity, 
clocks, and computation," Proc. of the IEEE, vol. 98, no. 
11, pp. 1828 -1846, nov. 2010. 

[2] Yik-Chung Wu, Q. Chaudhari, and E. Serpedin, "Clock 
synchronization of wireless sensor networks," IEEE Sig- 
nal Process. Mag., vol. 28, no. l,pp. 124-138, jan. 2011. 

[3] F. Gustafsson and F. Gunnarsson, "Mobile positioning 
using wireless networks: possibilities and fundamental 
limitations based on available wireless network measure- 
ments," IEEE Signal Process. Mag., vol. 22, no. 4, pp. 41 
-53,july2005. 

[4] R. T. Raj an and A.-J. van der Veen, "Joint ranging and 
clock synchronization for a wireless network," in Proc. 
4th IEEE International Workshop on Computational Ad- 
vances in Multi-Sensor Adaptive Processing ( CAMSAP), 
dec. 2011, pp. 297 -300. 

[5] S. P. Chepuri, R. Rajan, G. Leus, and A.-J. van der Veen, 
"Joint clock synchronization and ranging: Asymmetrical 
time-stamping and passive listening," IEEE Signal Pro- 
cess. Lett, vol. 20, no. 1, jan 2013. 

[6] Jun Zheng and Yik-Chung Wu, "Joint time synchroniza- 
tion and localization of an unknown node in wireless sen- 
sor networks," IEEE Trans. Signal Process., vol. 58, no. 
3, pp. 1309 -1320, march 2010. 

[7] Yiyin Wang, Xiaoli Ma, and G. Leus, "Robust time-based 
localization for asynchronous networks," IEEE Trans. 
Signal Process., vol. 59, no. 9, pp. 4397 -4410, sept. 
2011. 



